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^4—1 ■ Abstract 

o ■ 

[/J ' Phase behaviour of the Stockmayer fluid is studied with a method 

similar to the Monte-Carlo annealing scheme. We introduce a novel order 

Cu , parameter which is sensitive to the local co-orientation of the dipoles of 

S ' particles in the fluid. We exhibit a phase diagram based on the behaviour 

I I of the order parameter in the density region 0.1 < p* < 0.32. Speciflcally, 

' ^ . we observe and analyse a second order locally disordered fluid —>■ locally 

K ' oriented fl,uid phase transition. 

o- 
o 

^ 1 INTRODUCTION 

(N 

^ ' The study of dipolar spheres models is of significant importance as they can be 

considered as native models on ferrofluids [I] or dispersions of polar molecules, 
or as simple media with a property of self-organization, [5] . 



^_ ■ There are three models frequently studied in the context of dipolar flu- 

f— ^ I ids: the dipolar hard sphere (DHS) model, the dipolar soft sphere (DSS) model, 

psj ' and the Stockmayer (SM) model. The total potential energy of each of these 

models has the following form: 

_\^ _' with Udd being a dipole-dipole potential, Usr ~ a short range potential and 

summation being performed over all pairs of particles. 

The first term Udd is common for all three models; they differ with the respect 
to their short-range potentials: 
the DHS model has the hard sphere short range potential 

rjij ^ / 0, if ry > R, 

''^ \ oo, if nj < R. 
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the DSS model employs the soft sphere potential 

12 



the Stockmayer (SM) model is distinguished by its use of the Lennard- Jones 
short range potential 
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In these formulae, r^j stands for the distance between i-th and j-th particle. 
In the current paper we shall pay attention to the last model mentioned. 
The potential energy of the Stockmayer fluid of N particles is written as 
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Study of the DS models uses numerical methods since no theory for 
quantitative description has been developed. The Monte-Carlo and molecular- 
dynamical studies resulted in significant achievements in discovering the phase 
behaviour of these models. A large part of these results is summarized in a 
review by Weis and Levesquc, [3]. 

Thus, previous studies of this topic have shown that all three basic mod- 
els form an isotropic phase at high temperatures. At the same time, at high 
densities these models demonstrate low-temperature crystal phases, [HIS]. This 
difference of phases arouses interest in the phase transitions in these systems. 

Recent studies of phase transitions have revealed the presence of the 
isotropic liquid — > orientationally ordered liquid phase transition in dilute DS 
fluids. Specifically, at low densities the SM fluid exhibits the so-called string 
liquid phase, formed by chains of dipolar particles with dipoles co-orientated in 
each chain, see fig. [1] 

Despite all the similarities of the three models, the Stockmayer model 
demonstrates more complicated behaviour due to the Lennard- Jones potential 
it employs. Thus, the SM model demonstrates a phase transition gas — >■ 



isotropic liquid, [21IS1I3, that was not observed in the DHS/DSS models. 
Because of the interest in the magnetic properties of the Stockmayer model, 
a large part of its study has been performed with an external field applied. 
The introduction of an external field, however, originates a special direction in 
space which might have a bearing upon the behaviour of the system. 

Despite the great number of theoretical and simulation studies on the 
Stockmayer model, the full phase behaviour, including effects due to the 
entropy, has so far not been presented. We suggest a new technique of 
investigating based on measuring the local order in the DS models. It is 
worthwhile to note that the results of our calculations are in accordance with 
previous reports on the SM model. 

Study of the SM model might also shed light upon specific features of water. 
Each water molecule has a high dipole value that corresponds to dipoles in DS 
models. In the same time the Lennard- Jones potential is frequently used in 
molecular-dynamical models on water, e.g. TIP3P or TIP4P, and apparently 
approximates fairly the intermolecular interaction between its particles. The 
combination of these facts brings us to the idea of considering the Stockmayer 
fluid as an extreme model on water. 

Another fact that speaks in favour of this thesis is that the particles of 
both water and SM fluid contain the feature of self-organization into com- 
plex structures, such as chains (fig. [J) or clusters of particles, though the 
mechanisms of such organization are different. 







■■ '"V/, ' \ 1 ' \. *. V-''"" V 










i' "~\-~ ' ' -- \^ ^ " 


-'(',"■ i; '■^- ^!^- 






\ •- . ■■ ' r 1.--'^ ^ "*- 


-,' . ■- Z"--^'' 






















■..'*; ; -\ /V)^r^ ^^; r<\ 






















''','"'''' ''^v' . .-(''-'''" "■ -. \\/ '^'-y^^^ \ 


\{- ""' -■"i'"^~f ""' '\J \ I ""V'j: ^Z"" ""'""-'^.'^ 




















"^V- - - 1 " "^ \' ■ ), \ .^^"^ 





Figure 1: String phase in the Stockmayer fluid at low densities. Projections of 
dipoles on a plane arc shown. 



Let us consider a liquid of elongated molecules (e.g. rod or bar shaped), 
with its particles somehow ordered in space. If this order is stable, i.e. conserves 
for a long time period, this liquid is a liquid crystal (LC): it flows under shear 



stress and at the same time has the property of an anisotropic arrangement of 
molecules, which is typical for crystals. Depending on the character of such 
order, three types of LC phases are usually distinguished: nematic, cholesteric 
and smectic. 

In the nematic phase all molecules arc approximately directed along one 
line called the nematic axis. The microscopical structure of the other two 
is more complicated: these liquid crystals have a lamellar structure with 
molecules co-orientated along certain axis in each layer. The difference between 
the directions of the axes of neighbour layers is usually small. 



Considering a microscopical structure of liquid crystals, we would like 
to pay attention to an important property of local orientational order. In a 
perfect nematic phase every two molecules are almost co-directed, i.e. a global 
orientational order takes place. In a real liquid crystal, however, the property 
of co-orientation weakens with distance since the direction of the nematic axis 
may vary for different parts of nematic sample. 

In a substance with local orientational order the directions of two molecules at 
a short distance are correlated, and the correlation decreases with distance. A 
fluid with global orientational order is a locally orientated fluid, but not vice 
versa. 

During our calculations we have observed the similar effect of local ordering in 
low temperature phases of the Stockmayer fluid. However, in this case we dealt 
with the orientation of dipoles instead of such of molecules. The utilization of 
liquid crystal terminology produces similar terms for phases of the SM model. 
Thus, in paper |14| . phases of hard and soft spheres with orientational order 
are called ferroelectric nematic. The word ferroelectric underlines the dipolar 
meaning of the orientation of particles. 

Recent studies (as well as our calculations) confirm the presence of local 
orientational order in low-temperature phases of the Stockmayer fluid. The 
presence of local order, however, does not generally imply an existence of any 
order at global range, fig. [21 




Figure 2: A phase with local orientational order does not have global orienta- 
tional order. Projections of dipoles on a plane are shown. 

Since the presence of local orientational order is a distinctive feature of low- 
temperature phases, we shall make an attempt to describe phase behaviour of 
the Stockmayer fluid from this point of view. 

2 LOCAL ORDER PARAMETER 

According to the general theory of Landau - de Gennes, a second order phase 
transition can be described by means of order parameter. In our case this is 
a quantity that can put apart phases with local order and without it. Thus, 
the following discussion is performed in framework of the Landau - de Gennes 
theory. 

The main point of this work is our suggestion of a novel order parameter 
with the latter ability in the following form: 
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Any pair of co-directional dipoles in SM fluid increases the order parameter 
just as any opposite-directed pair decreases it. Moreover, this influence weakens 
rapidly with the distance increasing between the elements in such pairs. As 
a result, the value of G depends basically on the close-ranged orientation of 
dipoles in the SM fluid. In this manner the suggested order parameter detects 
a local order of the model and neglects any long-range disorder. 



It is worth noting that this order parameter is not the only one that 
takes into account the local order of particles. The one we use is simple to 
count and therefore rather practical. 



It is worth noting, that in spite of the fact that the presence of local 
orientational order in the Stockmayer model has been reported in previous 
studies, we could not find any quantitative description of this property crucial 
for the analysis of the phase behaviour that we performed. 

3 SIMULATION METHOD 

We performed a Monte Carlo (MC) calculation of a cubic box with Stockmayer 

particles in an NVT-ensemble for an interval of temperatures with fixed these 

basic parameters: 

Lennard- Jones depth e is 0.25 for all particles, 

dipoles value D is 5 for all particles, 

particle radii - R - are 1, 

particles number N is 2000 or 5000 (all calculations were made for both 

parameters) . 

Since we have to perform a Monte-Carlo simulation for a large range of 
temperatures, there is a need for accelerating our calculations. To overcome 
this difficulty we apply a method that is similar to the MC annealing scheme, 
[13] . The algorithm starts from a random particles distribution at high tem- 
perature (i.e. in a gas phase). After about 1-^2- 10* MC steps the calculation 
stops and the current system conformation becomes a start condition for the 
new computation with a slightly decreased temperature. This technique allows 
us to reduce a number of MC steps per temperature and thus significantly 
reduce calculation time. 

It is also worth noting that the annealing method is used for finding the 
global minimum of a function of several parameters. Thus, we may expect to 
see the system we simulate in the lowest energy state at zero temperature. 
Such behaviour is expected from a real system as well. This fact supports this 
technique. The process we simulate can thus be corresponded to slow cooling 
of a real system. 

4 RESULTS 

We shall utilize the notation used in paper [3] for the reduced density p* = ^^^ 
with V - cell volume. A'' - number of particles and R - particle radius. 
The temperature scale T* we shall use is based on the following estimates. If 
colloid particles have radii of order ~ lO^A = 10~^cm and the magnitude of its 
dipoles is ^ 10 Debye, we can thus roughly estimate the dipolc-dipole interaction 
energy as 



L>2 100 
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Using this value as an estimate of the mean energy of a single particle in the 
fluid {kT* - Udd) we obtain for T=300K 



kT 
kT* 



1,38-10 



-16 



300 



10 



-13 



0,41 



Thus, 0, 41/cT* corresponds to the temperature of 300K. Despite the roughness 
of these estimates, the order of its magnitudes indicates that the results we 
obtain are physically adequate. 

Let us now consider a sample of the typical behaviour of the local order 
parameter G and the potential energy of the system. The figures [3] and |4] 
contain temperature dependencies of these quantities. 
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Figure 3: Temperature dependence of order parameter G. Reduced density p* = 
0, 16, number of particles N=5000. 

We shall pay our attention to the following reduced temperature regions in 
the graph [31 

A : 0.01 - 0.42 : corresponds to a phase with a local orientational or- 
der (dipoles of the SM particles are co-orientated at close range) ; 

B : 0.42 - 0.65 : corresponds to a zone of phase transition. The be- 
haviour of the order parameter is similar to that of magnetization in the Ising 
model and indicates a second order transition (in Erenfest's classification); 

C : 0.65 - 1.40 : is the range of the locally (and globally) orientation- 
ally disordered phase. 



It is necessary to make a correction in the classification above. There is 
one more interesting point c' at kT* « 0,9. This point is hard to see on the 
order parameter graph [3] because of the large scale, but one can notice it on the 



energy graph U) At point c' both graphs are bent, or speaking mathematically, 
have a discontinuity of first derivative (not in a precise sense however). Thus, 
this point corresponds to the second order gas -^ liquid phase transition. 
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Reduced temperature, kT* 

Figure 4: Temperature dependence of potential energy. Reduced density p* = 
0, 16, number of particles N=5000. 

From this point of view, the B region (see above) corresponds to the 
disordered liquid — > liquid with local orientational order phase transition. 

There are two additional points in the fig. [3] that arouse interest: kT* k, 0.46 
and kT* k, 0.55. These points correspond to the local maxima of the first 
derivative of G. We shall recall them in the further text. 



We would like to underline the fact that the potential energy behaviour 
does not provide phase information as clearly as the order parameter does. 
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Figure 5: Phase diagram of the Stockmayer fluid based on the local order pa- 
rameter. Solid lines a, b, d correspond to conformational transitions. Dashed 
line c designates local maxima of first derivative of the local order parameter 
G. 

Let us now consider phase diagrams based on the behaviour of order 
parameter, fig. [SJ [SI The graph on fig. [3] is a vertical section {kT* — const) 
of the phase diagrams [5] and HI The regions A, B, C of the graph [3] now 
correspond to large zones of phases, bounded by the lines d and b (fig. [5]). The 
point c' correspond to the a line. Finally, local maxima of order parameter G 
correspond to dashed c line. 



Thus, we have 3 basic lines (solid lines a, b, d on fig. [5]), that perform 

the division of the phase diagram on 4 parts: 

the zone above the d line corresponds to a phase with local orientational order; 

in the region between the b and d lines there is a transition zone (dipoles 

become locally co-orientated); 

in the region between the b and a lines the SM fluid has an isotropic liquid 

phase; 

the region below the a line corresponds to an isotropic gas phase. 
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Figure 6: Phase diagram of the Stockmayer hquid based on the local order 
parameter. Four basic phases are distinguished. 

Our calculation shows that conformational organization of the Stockmayer 
particles depend both on temperature and the reduced density p* . Here we 
perform the analysis of conformations of particles in phases at different temper- 
atures and densities. To that end we shall consider states of the system on the 
last step of Monte-Carlo algorithm, calculated for different values of kT* and p* . 

In the lowest zone (below the a line) on figure [5] the fluid exhibits a uni- 
form isotropic gas phase for all considered values of density. The region 
between the a and the b lines demonstrates an isotropic phase with density 
fluctuations. Based on the presence of clusters we regard this phase as liquid. 
The decrease in temperature above the b line results in the appearance of 
a weak local order of dipoles. At this point dipolar vortexes and regions of 
co-orientation of dipoles appear in dense parts of medium. The cooling of the 
system from the b to the d line contains a cascade of two-step reconstructions 
of the fluid: (1) particle displacement and formation of a new topological 
structure followed by (2) dipoles re-organization induced by new positions of 
particles. 

The second part of each of these reconstructions corresponds to a local maxi- 
mum of first derivative of local order parameter G, points that we have shown 
as a dashed line on the fig. [5l The double dashed line in the high-concentration 
region stands for two reconstructions we managed to distinguish in this interval 
of the reduced densities. 

We feel that this cascade is an interesting effect due to entropy, and it needs to 
be examined in more detail. 

Thus, during our calculation we observed a reconstruction that is reminiscent 
of the BKT transition in XY-model: a phase with vortexes transformed to a 
phase with linearly co-orientated parts in a very short range of temperature 
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(about OmkT*). 

The Stockmaycr fluid phases at low temperatures (above the d hne) are 
locally orientationaUy ordered. In this region of the phase diagram we could 
distinguish three basic types of low-temperature conformations: liquid of clots, 
'network' phase and phase of opposite-directed nematic strands, see fig. [71 
It is worthwhile to note, that the existence of the latter phase has not been 
reported before, in spite of the vast number of publications on the topic. We 
were unable to determine the exact concentration boundary between these 
low-temperature conformations; these transitions could be gradual. This region 
of phase diagram is to be investigated in more detaiL 



= 0.1 



(o) 







Figure 7: Schematic images of three low-temperature phases. Arrows designate 
local directions of dipoles. The first picture corresponds to a liquid of globules 
with dipoles forming vortexes. The second picture corresponds to a complex net- 
work structure. The third picture represents opposite-directed nematic strands. 
Periodic boundary conditions are implied. 

5 CONCLUSIONS 

Despite the large number of theoretical and computational studies, we are 
still far from understanding the phase behaviour of the Stockmaycr fluid. 
The complexity of effects due to the entropy in the model and flaws of both 
theoretical and computational methods of research demand developing novel 
techniques in the study of the Stockmaycr fluid. 

In the current paper we suggest such a technique, which is based on a notion 
of local orientational order. The property of local orientational order is an 
important characteristic of low-temperature phases of the Stockmaycr fluid. 
We have suggested a novel order parameter G as a measure of such order. The 
utilization of this parameter allowed us to perform an analysis of the phase 
behaviour of the Stockmayer fluid from a new point of view. 
Thus, using our technique resulted in the construction of the phase diagram 
(fig. [5]). We also managed to describe the isotropic liquid -^ locally ordered 
liquid phase transition and find a low-temperature phase, that composed of 
opposite directed parallel nematic layers (fig. [71 right). The full phase diagram 
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(for < p* < 1) and analysis of conformations for different values of kT and p* 
shall be presented in our further publications. 

The local order we introduced provides a new point of view on systems without 
global order. It might be useful to discover hard-to-find inner symmetries in 
other soft matter systems. Specifically, we expect this method to be useful in 
models related to liquid crystals or water. It might also provide new results in 
self-assembly models. 

During our research we have also used a Monte-Carlo annealing tech- 
nique that turned out to be a very effective numerical instrument for studying 
soft matter systems. 

This work was supported by RFBR, grants #ll-02-01462-a, #10-01-00748 
and by the Government Grant of the Russian Federation for support of 
research projects implemented by leading scientists, Lomonosov Moscow State 
University under the agreement No. 11.G34.31.0054. 
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